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Hoofdstuk I. Inleiding 


Deze handleiding is bedoeld voor Algol-kenners die in 
Fortran willen gaan programmeren. Omdat dus niet ingegaan 
hoeft te worden op de meest elementaire begrippen van het pro- 
grammeren is het mogelijk wat meer te behandelen dan in de 


FORTRAN-cursus(ACCU reeks nr. 2). Bovendien wordt speciaal de 


nadruk gelegd op punten, die in Fortran zodanig verschillen 


van Algol, dat ze veelvuldig aanleiding tot moeilijkheden 

geven. De nummers tussen haakjes bij de paragraaftitels ver- 
wijzen naar de betreffende paragrafen van de cursus Algol 60 
(A. v.d. Sluis en C.A.Ch. Görts; Cursus Algol 60). | | 

Er bestaan verschillende dialecten van Fortran, waaruit 
een ASA-commissie een Standaard Fortran heeft samengesteld. 
We zullen ons hier beperken tot dit Standaard Fortran (ook 
wel Fortran IV genoemd) en slechts enkele extra faciliteiten 
vermelden die op de Cyber mogelijk zijn. Als men echter deze 
faciliteiten gebruikt, moet men er rekening mee houden dat 
het programma waarschijnlijk niet zonder meer op een machine 
van ander fabrikaat verwerkt kan worden. | 
De Cyber kent twee Fortran-compilers FIN en RUN (zie ACCU MAP 
bulletin ff, 3). 

In deze syllabus wordt niet ingegaan op mogelijkheden 
van Fortran, die over het algemeen slechts bij geavanceerde 
programmeerproblemen optreden. Voor deze gevallen en ook voor 
specifieke faciliteiten op de Cyber raadplege men: 

CYBER 70 FORTRAN EXTENDED REFERENCE MANUAL (FTN) 
CYBER 70 FORTRAN REFERENCE MANUAL (RUN). 


Daar men een Fortranprogramma in het algemeen op pons=- 
kaarten aanbiedt, willen we hier iets zeggen over de kaart- 


indeling. Een ponskaart heeft 80 kolommen. Hiervan zijn alleen 
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de eerste 72 van belang voor een programma. De kolommen 73 

t/m 80 kunnen gebruikt worden om de kaarten te nummeren. 

(In datakaarten kunnen alle 80 kolommen gebruikt worden.) 
Er geldt: 


als in kolom 1 een C staat wordt de rest van die kaart 

als comment beschouwd. | 

de kolommen 1 t/m 5 (in kolom 1 geen C) kunnen een 
statement-nummer (gelijk aan label bij Algol) bevatten. 
Dit is een getal van 1 t/m 5 cijfers. 

kolom 6 is de vervolgkolom. Als een statement, declaratie 
of specificatie te lang is voor één kaart mag hij vervolgd 
worden op de volgende kaart. Men moet dan een willekeurig 
karakter , ongelijk aan een blank of nul, in kolom 6 van de 
vervolgkaart(en) plaatsen. In de eerste kaart van een 
statement, declaratie of specificatie moet kolom 6 een 
blank of een nul bevatten. Na een commentkaart mag geen 
vervolgkaart komen. Voor meerdere regels comment moet 

dus iedere kaart in kolom 1 een C bevatten. 

de kolommen 7 t/m 72 bevatten de programmatekst. Blanks 
hebben geen invloed op het programma en kunnen willekeurig 
tussengevoegd worden om de leesbaarheid van de print-out 
te verbeteren. | 

op iedere kaart staat niet meer dan één statement.(Op de 
Cyber is het echter wel mogelijk meerdere statements per 
kaart te ponsen. Ze moeten dan gescheiden worden door 


een $ dollar-teken). 
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I.3 Karakterset 


Fortran kent o.a. geen kleine letters en ook de onder- 
streping bestaat niet. De karakterset bestaat uit: 
a) de letters: 
ALG BEE GH a ute e à «eht Á 
b) de cijfers: 
024456 sd 9 
ec) de speciale karakters: 
blank (spatie) + - / +.) ,(!' $ 


Hieronder volgt van een aantal Algol-symbolen het equivalent 


in Fortran: 


ALGOL (op band) FORTRAN (op kaart) 
n Bi. 

es NE 

< ek 

< ‚LE. 

> Gls 

> ‚GE. 

Á .AND. 

V OR: 

| ‚NOT. 

X # 

A kk 

[ ( 

] ) 

zo (laag tientje) E 

8 einde kaart of $ 
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Een identifier bestaat uit 1 t/m 6 letters en cijfers 


en moet beginnen met een letter. 


Een statementlabel is een getal van maximaal 5 cijfers. 


Het statementnummer moet in de kolommen 1 t/m 5 geponst 


worden. 


Ga) 


Gb) 
Ge) 
Gd) 
Ge) 


GÉf) 


Gg) 


Getallen (20) 


Een getal heeft een van de vormen: 


(ddd is een rwijtje cijfers; minstens Één) 


ddd | (integer) 
ddd. | 

Boes (real) 
ddd.ddd 


Een van de mogelijkheden Ga, Gb, Ge of 
Gd gevolgd door: E ddd of E + ddd of 
E - ddd 
Een van de mogelijkheden Ga, Gb, Ge of (double precision) 
Gd gevolgd door: D ddd of D + ddd of 
D — ddd 
(E en D komen overeen met laag 1o=tje bij 
Algol) 
(real,real) (complex) 
Zo wordt b.v. het complexe getal 1+3.4 i 
aangegeven door (1.,3.4). Het is aleen 
mogelijk om van twee reële getallen op deze 
wijze een complex getal te maken. Indien men 
twee reële variabelen tot een complexe groot- 
heid wil samenvoegen moet men de standaard- 
Luhette CMPLMÁ saasrrad Bebruiker: 


Opmerking: Als A b.v. de waarde 10-5 moet krijgen, 


hek Ein 


moet dat gebeuren met A=1.E-5 en niet met 
A=E-5 want dit laatste betekent de variabele 


E min 5. 


Een string kan geschreven worden met behulp van de 


H-descriptor (zie VIII.4). Op de Cyber is ook toegestaan 
USS... S% 


hierbij is s een willekeurig karakter niet gelijk aan %#. 
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Hoofdstuk III. Variabelen en hun declaraties 


Fortran kent de z.g. impleciete declaratie, d.w.z. als 
niet expliciet'anders gedeclareerd wordt zijn variabelen, 
die met een I, J, K, L, M of N beginnen, van het type INTEGER, 
terwijl variabelen, die met een andere letter beginnen van 
het type REAL zijn. In Fortran hoeven variabelen dus niet 
altijd gedeclareerd te worden. Als men hiervan af wil wijken 
is er wel een expliciete declaratie nodig. Deze moet net als 


bij Algol aan het begin van de betreffende programma-eenheid 
staan. | 


Een declaratie voor een eenvoudige variabele heeft de 
vorm: 


INTEGER Wbs 
REAL Unen 
DOUBLE ERECISLON: WE ssasstl 
COMPLEX ON enn 
LOGICAL VIseensdd * (Im Algol: Bóolean) 


(VI betekent Variabele Identifier) 


Een DOUBLE PRECISION grootheid beslaat twee geheugenwoorden 
en heeft op de Cyber een nauwkeurigheid van ongeveer 28 deci- 


male cijfers (gewone REAL 14 cijfers). 


III.3 Array declaratie (23) 


Een array declaratie is van de vorm 


AD AI(n),...,AI(n) 


waarbij 


AD is 1) een van de type aanduidingen van II.2. 
2) DIMENSION 
3) COMMON (voor betekenis van COMMON zie VI.6) 


en 


Alin) Ìe 1) Al(n,) 
a) Aln, ‚n,) 
3) Aln, ‚n,,n3) 


Hierbij zijn n,» nj en nz gehele getallen (of 
integer variabelen, zie III.5). 


(AI betekent Array Identifier) 


De indices lopen van 1 t/m n 1 t/m n, resp. 1 EAM Mio» 


1 3 


Omdat de indices dus altijd bij 1 beginnen wordt geen 
_benedengrens aangegeven. Er zijn maximaal slechts 3 indices 
mogelijk. 

Het type van de indices moet INTEGER zijn. 


Voorbeeld: 


DIMENSION A(2,3), B(6), C(6) 
Nu is A een reële matrix, waarbij de eerste index de waarde 


1 en 2 aan kan nemen en de tweede index van 1 t/m 3 loopt. 


N.B. Let op het gebruik van ronde haken; vierkante haken 


bestaan niet in Fortran. 


Opmerking: Achter iedere array-identifier moet de bovengrens 
expliciet worden aangegeven. We kunnen in het bovenstaande 
voorbeeld de array declaratie voor B en C dus niet samen- 
trekken tot B‚C(6). | 


De array-elementen krijgen in het geheugen op elkaär 
volgende plaatsen toegewezen. De volgorde wordt bepaald door 
de regel dat de eerste index het eerste toeneemt en de laatste 


index het laatste. 
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B.v. het array A met 
DIMENSION A(3,2,3) 
komt in het geheugen in de volgorde 


AC1,1,1) A(2,1,1) A(3,1,1) A(1,2,1) A(2,2,1) A(3,2,1)> 
sAllsld) ACL ABL) Ald) ACH AIDE 
sACTsL, 3) Meda) AAA, Ald) ALLA ACID) 


TIt.S Variabele dimensië 


Als minstens n; in een array-toewijzing een integer 
variabele is, heet het array een "variabele array!" en de 
variabele wordt een "variabele dimensie" genoemd. Een variabele 
array mag slechts in een subprogramma optreden en de naam van 
het array en de variabele dimensie moeten tot de formele para- 
meters van dit subprogramma behoren. In het programma, waar 
dit subprogramma aangeroepen wordt, moet aan het array,dat 
als actuele parameter in de aanroep voorkomt ,een array=- 
toewijzing met vaste grenzen zijn gedaan, die niet kleiner is 
dan de toewijzing die volgt uit de variabele dimensie. Dit is 


dus geen "echte" variabele dimensie zoals Algol die kent. 


Een statement waardoor men invloed kan uitoefenen op de 
geheugentoewijzing is de EQUIVALENCE-statement. Hierbij wordt 
aan de in de EQUIVALENCE-lijst verbonden grootheden dezelfde 


geheugenplaats(en) toegewezen. Zo wordt met 


DIMENSION E(5) | 
EQUIVALENCE (A,B,C),(D,E(3)) 


aan de grootheden A, B en C éénzelfde geheugenplaats toegewezen. 


Evenzo wordt aan D en E(3) één (andere) geheugenplaats toege- 
kend. | 


Nog een voorbeeld 


COMPLEX C 
DIMENSION A(2) 
EOUIVALENCE (CA) 


Nu bevat A(1) het reële deel van C en A(2) het imaginaire 
deel. | OO 
De EQUIVALENCE-statement wordt meestal gebruikt om 


geheugenruimte te besparen. 


Opmerking: Als in een EQUIVALENCE alleen de array=-naam zonder 
indices voorkomt, komt dit overeen met het eerste array- 
element. In alle gevallen blijft echter de volgorde in het 
geheugen, zoals besproken in III.4, gehandhaafd. 
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Hoofdstuk IV. Expressies 


In Fortran vertonen arithmetische expressies in het 
algemeen dezelfde gedaante als in Algol. Er worden echter 
voor vermenigvuldigen en machtsverheffen andere karakters 


gebruikt: 


vermenigvuldigen « (Algol op band x) 


machtsverheffen *k (Algol op band 4) 


Bovendien bestaat er zeen apart karakter voor de gehele 
deling, maar de gewone deling van twee integer grootheden 
heeft hetzelfde effect. Hier moet men als vroeger Algol- 
gebruiker bijzonder oppassen. Als b.v. A, B en C reêle 


grootheden zijn dan heeft 
Â = 3/2 
s daf 2 
Sal 2 


tot gevolg dat A=1. en B=C=1.5 
Men kan van een integer grootheid een reële grootheid maken 
met behulp van de standaardfunctie FLOAT. Zo heeft 


A = FLOAT(3) / FLOAT(2) 


wel de waarde 1.5 
Van een reele grootheid maakt men een integer met IFIX. 

De expressie voor B is een z.g. “mixed mode!" expressie. 
Volgens de ASA standaard van 1966 is alleen ‘mixed mode!" van 
reêle en dubbele precisie grootheden toegestaan voor de 
operatoren +, =, en /. De huidige compilers hebben echter 
dienaangaande over het algemeen geen beperkingen meer. 

Het type van het resultaat van een bewerking met 2 


operanden kan in het volgende schema samengevat worden 


operand 2 


de 
f 1 r dp G 

O 1 if  d dp c 

Pp 

e r r r dp a 

r 

a & dp dp dp dp ë 

n 

d c c c c c 

1 


Hierbij geldt: 


= integer 
r = real 
dp = dubbele precisie 
c = complex 


Het assigneren van grootheden van verschillend type aan elkaar 
gaat over het algemeen hetzelfde als in Algol. Een belangrijke 
uitzondering is echter het feit dat er bij toekenning van een 
reële grootheid aan een integertruncatie en geen afronding 


optreedt. 


Voorbeeld: 


Bij machtsverheffen geldt het volgende schema: 
Ax B 


rader de ENE dine nee Ed itn add ln TE Si EE TEE EE PE DE Rt EP EEE PPE 


kde a ke kn 








waarbij — betekent niet toegestaan. 


B.v. A=10##(-3) geeft Az=0 ! 


A VRIEND Ee EN 


De relaties worden geschreven met: 


| ‚GT. voor > 
| ‚GE. voor > 
| EO. voor = 
‚NE. voor $ 
„bl. voor < 
SDE hf voor < 


De logische operatoren zijn: 


ln an Ee et Ten ne a en a ie En ed lln 


| ‚NOT. in plaats van | 
| | „AND. in plaats van A 


.OR. in plaats van V 


Opmerking: = en > bestaan niet in Fortran. 
De logische waarden zijn .TRUE.; en .FALSE. 
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Hoofdstuk V Statements 


Een belangrijk verschil met Algol bestaat in het feit 


dat Fortran geen compound statement kent. 


De assigned statement heeft de vorm 


AE 
ln jn 
Opmerkingen: i) het =teken komt overeen met := in Algol 
ij) er mag volgens de Standaard maar één =teken 
per statement voorkomen. ko 
Op de Cyber zijn echter wel meerdere =tekens 
per statement toegestaan. Hierbij wordt de 
assigned statement van rechts naar links 


afgewerkt. | 
B.v. Az=lIl=B=3.4h wordt Bz3.b I=B AzlI 


De DATA statement dient om grootheden een beginwaarde toe 
te kennen. Dit geschiedt gedurende de compilatie. 
Zo wordt met | 


DATA Tad / dt k vAsBee Jl 2wbrsbs / 


aan I een beginwaarde 1 
J di 2 
B hi 5 
ook ; B 4 6 zit on 
en u mn 6 gegeven. 


Opmerking: Grootheden in COMMON kunnen niet op deze wijze van 


een beginwaarde voorzien worden. 
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V.3 De GO TO statement (30) 


De GO TO statement heeft de vorm 
GO TO n 


waarbij n een statementnummer 1s. 


Een computed GO TO statement heeft de vorm 
GO TO (n‚ ‚nos ea Ndsi 


met n. statementnummer en 


1 een integer variabele. 


Als 1 een waarde j met 1<j<m heeft, wordt het programma 
voortgezet met de statement met statementnummer n.. In de 
andere gevallen is het statement niet gedefiniëerd. De 


computed GO TO statement komt overeen met de Switch in Algol. 


De dummy statement luidt: 
CONTINUE 


Hij kan gebruikt worden voor het aanbrengen van statement- 


nummers , meestal aan het einde van een DO-loop. 


De conditionele statement (logische IF) heeft de vorm 
IF (BE) 8 


S 1s een uitvoerbare statement; het mag echter geen DO of 
logische IF zijn. 


Als BE true is wordt S uitgevoerd anders niet. 


Opmerking: Dit is de enige vorm van de logische IF statement 
in het Standaard Fortran. 
Op de Cyber is cok de vorm IF(BE)n, ‚n, mogelijk. Als BE true 


is, wordt naar label n, gesprongen, als BE false is naar n- 


1 


V.7 De arithmetische IF statement 


De arithmetische IF statement heeft de vorm 


IF (AE) n,snosN3 


met AE een arithmetische expressie van type integer, real of 
double precision en n‚.nj en n statementnummers. | | 


| 3 
Als AE < O0 wordt het programma voortgezet bij statementnummer 


Dj» als AE = 0 wordt het programma voortgezet bij statement- 
nummer n, en als AE > 0 wordt het programma vervolgd met 
statementnummer nn. 


9 


V.8 De DO statement (33) 


Een DO statement heeft de vorm 


DO n 1=M, >M, ‚Mg 


waarbij de DO-loop range zich uitstrekt tot en met het statement 
met het statementnummer n | 
1 is een integer variabele (DO-loop index) 


Mm, sm en m 


1m) zijn positieve integer getallen of integer 


3 
variabelen met positieve waarden. 


In ls de beginwaarde van 1 


1 
Mm, ls de eindwaarde van 1 
mg is de stapgrootte | 
Als m, gelijk is aan 1 mogen we verkort schrijven 


3 


DO n LEM, oM 


De DO statement komt overeen met de for statement bij Algol. 
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Voorbeeld: 


Bepaling van het grootste element in een array van 8 elementen 


DIMENSION A(8) 


DO 4 1=2,8 

IF (A(I)-G) t‚t,3 
3 G=A(I) 
hb CONTINUE 


Opmerkingen: 


1) 


2) 


3) 


5, 


à 


6) 


Een statement met label n moet na de DO volgen. Dit eind- 
statement mag geen statement van het sprong type zijn 
(geen GO TO, arithmetische IF, DO, STOP, RETURN of een 
logische IF met een van deze statements). 


Als de DO-loop verlaten wordt doordat i>m, is geworden, 


is 1 niet gedfiniëerd. Daarentegen kaa i zijn laatst 
toegevoegde waarde als de DO-loop door een GO TO of een 
arithmetische IF wordt verlaten. (Identiek met Algol). 
Verandering van de DO index 1 of de DO-parameters mm 


2 


en m, door statements binnen de DO-loop is niet geoor- 


3 
loofd, | 
Een sprong van buiten een DO-loop naar binnen is niet 
toegestaan. n OO 

Er zijn nesten van DO-loops mogelijk. (Ze mogen met het- 
zelfde label eindigen). 

Als m >, is wordt de DO-loop op de Cyber toch eenmaal 


dE 
doorlopen. 


Voor ‘declareren! van functies en subroutines zie 


Hoofdstuk VI. 


Het aanroepen van een subroutine gaat met 


CALL s (a, ; „a) 
of CALL s 


waarbij s de naam is van een subroutine en a; de actuele 
parameters zijn. 

Het aanroepen van functies geschiedt op dezelfde wijze | 
als in Algol. Voor de functienaam geldt weer dezelfde impliciete 
type aanduiding als voor identifiers. Een functie moet altijd | 


minstens één parameter hebben. 


V.10 De END statement 


De vorm van de END statement is 
END 


De END statement betekent voor de compiler het einde van een 
programma-eenheid. De END statement mag geen statementnummer 
hebben. Iedere programma-eenheid moet worden afgesloten door 


een END statement. 


V.11i De STOP statement 


De vorm van de STOP statement is 
STOP 


Uitvoering van dit statement heeft de beêindiging van het 
programma tot gevolg. Hij mag overal in het programma voor=- 


komen. 
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Hoofdstuk VI. Statementfuncties, functiesubprogramma's en 


subroutines 


De verbindingen tussen de subprogramma's en het hoofdpro- 
gramma zijn in Fortran heel anders dan bij Algol. Subpro- 
gramma's zijn afzonderlijke programma's (programma-eenheden) 
en hun enige verbinding met het programma waarin ze aange- 
roepen worden bestaat in de meegegeven actuele parameters 

en de verbinding die door de variabelen in COMMON (zie VI.6) 
tot stand wordt gebracht. De andere variabelen zijn locaal. 


Recursiviteit,ook in verkapte vorm, is niet toegestaan. 


VI.1 De Statementfunctie 


Een statementfunctie wordt aan het begin van een pro- 
gramma-eenheid, direct na de declaraties, door een enkele 
omschrijving, die de vorm van een arithmetische of logische 
toekenning heeft, gedefiniëerd. Een statementfunctie is een 
functie-procedure waarin de arithmetische expressie in één 
statement kan worden opgeschreven. De aldus gedefiniëerde 
statementfunctie is een intern onderprogramma van de gegeven 
programma-eenheid en mag alleen in deze programma-eenheid 
gebruikt worden. De naam van een statemenfunctie mag niet in 
een DIMENSION, EQUIVALENCE, COMMON of EXTERNAL beschrijving 
optreden en mag in de programma-eenheid waarin ze voorkomt 


niet als naam van een variabele of array gebruikt worden. 


Voorbeeld: 


Als een programma een aantal keren de wortels uit vierkants- 


vergelijkingen berekend moet worden, kunnen we schrijven: 


ROOT1(A,B,C)=(-B+SQRT(B&4&2-URARC) )/(24A) 
ROOT2(A,B,C)=(-B-SQRT(Be&2-hwARC))/(24A) 


. e e e Ee 


X1=ROOT1(F,G,H) 
X2=ROOT2(F,G,H) 





rd 
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Functiesubprogramma!'s (Algol: functie procedure) 


Een functiesubprogramma is van de vorm 
E FUNCTION fla,s...sa) 


programma-body (Hierin moeten voorkomen 
f=... en RETURN) 


END 


type aanduiding INTEGER, REAL, DOUBLE PRECISION, COMPLEX, 
LOGICAL; (Deze mag worden weggelaten als de naam f in 
overeenstemming is met de impliciete typt aanduiding) 
symbolische naam van de functie (met dezelfde eisen als 
voor identifier naam) 

formele parameters; namen van variabelen, namen van 
arrays of namen van functies. De specificatie van de 
formele parameters geschiedt op dezelfde wijze als in 

een hoofdprogramma. Een array moet met grenzen aangegeven 


worden; dit is dus afwijkend van Algol. 


Er moet gelden: 


1 


2) 


3) 


U) 


oi, 


6) 


Bij het doorlopen van het functiesubprogramma moet aan 
de naam van de functie, geschreven als variabele zonder 
parameters, minstens eenmaal een waarde toegekend worden. 
De naam van de functie mag in het functiesubprogramma in 
geen declaratie statement voorkomen. 

De namen van de formele parameters mogen in het functie- 
subprogramma niet in EQUIVALENCE of COMMON voorkomen. 
Een functiesubprogramma mag zichzelf in zijn eigen 
definitiegebied niet direct of indirect aanroepen (geen 
recursiviteit). 

In het functiesubprogramma moet minstens één RETURN 
statement voorkomen. 


Er moet minstens één formele parameter zijn. 


Een functiesubprogramma wordt aangeroepen doordat men ze als 


primary in een arithmetische of logische expressie gebruikt 


voorzien van een geëigende lijst van actuele parameters 


(overeenstemmend met de formele naar aantal, type en volgorde). 


De actuele parameters mogen zijn: 


1) constanten 

2) namen van variabelen 

3) namen van arrays of arrayelementen 

kb) een willekeurige expressie 

5) namen van andere functiesubprogramma's en sub- 


routines. 
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(Deze namen moeten dan ook in een EXTERNAL declaratie 


voorkomen (zie VI.5)). 


In appendix A geven we een lijst van standaard=-functies die 





zonder definitie gebruikt mogen worden. 





VI.3 Subroutines 


Een subroutine heeft de vorm 


SUBROUT INE TC, 

of SUBROUTINE s 
programma=body met minstens Één keer een RETURN 
END | 


symbolische naam 





as: formele parameters 


Daar een subroutine sterke overeenstemming vertoont met een 





functiesubprogramma gelden met uitzondering van 1) en 6) 
dezelfde eisen. 
Het aanroepen van een subroutine gebeurt door een CALL 


statement. 


CALL s(b,»---sb_) 
Gt CALL s 


ds ikk ie nn Se nk elan et da ann meere d merken TE AE 


waarbìj de actuele parameters (indien aanwezig) weer in 
aantal, volgorde en type moeten overeenstemmen met de formele 


parameters. 
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a) 


b) 


DA 


Voorbeelden: 


Subroutine ter bepaling van de G.G.D. van twee integers 


SUBROUTINE GGD(I,J,K) 

I=IABS(I) 

J=IABS(J) 

TECO Zed 
1 L=I-(I/J)#J 

I=J 

J=L 

TRC 4e0d 
2 K=I 
RETURN 

END 


Om de G.G.D. van b.v.v 21 en 56 te bepalen moeten we deze 


subroutine aanroepen met 
CALL GGD(21,56,M) 


M krijgt dan de waarde 7. 

Omdat bij het bepalen van de G.G.D. slechts één waarde 
teruggegeven moet worden aan het aanroepende programma 
kan dit natuurlijk ook in de vorm van een functiesub- 


programma gegoten worden: 


INTEGER FUNCTION GGD(I,J) 
I=IABS(I) 
ze 
2 GGD=I 
RETURN 
END 


Met nu als aanroep 


M=GGD( 21,56) 


Als tweede voorbeeld geven we de matrix-vermenigvuldiging 
van vierkante matrices. We gebruiken hierbij de z.g. | 


variabele dimensie. 











DIMENSION E(5,5),F(5,5),G(5,5) 


hoofdprogramma 


SUBROUTINE MATMUL(A,B,C‚‚N) | 
DIMENSION A(N‚,N),B(N,N) ,C(N‚N) 
DO 2 I=1N 
DO 2 J=1N 
S=0 
DO 1 Kz=1N 
1 Sz=S + A(I,K)#B(K,J) 
% GETS 
RETURN 
END 


subroutine 


VI.b De RETURN statement 


De REIURN statement bewerkstelligt de terugkeer naar het 
aanroepende programma. Hij mag overal in een subprogramma 


voorkomen. 


VI.5 De EXTERNAL statement 


@ Kd Kd V 
zee RE 


v;: naam van een subprogramma. 


EXTERNAL v, ‚v 


De EXTERNAL statement moet gebruikt worden wanneer een 
subprogramma als actuele parameter in de aanroep van een ander 
subprogramma optreedt. De EXTERNAL statement moet staan aan 
het begin van de programma-eenheid waarin deze aanroep plaats- 


vindt. 


Voorbeeld: 


sen 


EXTERNAL POLY 


RETURN 
END 


FUNCTION POLY(X,Y,Z) 


POLLS ada 
RETURN 
END 


hoofdprogramma 


subroutine 


functiesubprogramma 


Wat de verbindingen tussen de verschillende programma- 


eenheden betreft zijn er grote verschillen tussen Algol en 


Fortran. Fortran kent geen blokstructuren zoals Algol. Elke 


programmaeenheid (hoofdprogramma, functiesubprogramma of 


subroutine) wordt in Fortran als een afzonderlijk programma 


gezien, dat alleen verbindingen heeft met een andere pro- 


gramma-eenheid door de actuele parameters en door de groot- 


heden die 


in COMMON staan. 


Een COMMON statement heeft de vorm 


COMMON Vis. esVE 


De COMMON statement dient-om in een programma-eenheid 


grootheden uit andere programma-eerheden (eventueel zelfs 


onder een andere naam) ter beschikking te krijgen zonder dat 


deze grootheden in de parameterlijst voorkomen. Ze moet aan 
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het begin van de betreffende programma-eenheden staan. 
Voorbeeld: 


In hoofdprogramma 
COMMON A,B‚C,‚D(2) 
(hierbij is D(2) ook een array declaratie, zie II.3) 


en in subprogramma 
COMMON P,A,C,R(2) 


Dan wordt dezelfde geheugenplaats toegewezen aan 
A uit het hoofdprogramma en P uit het subprogramma 


en evenzo aan 


B uit het hoofdprogramma en A uit het subprogramma 
C LÂ Á LÌ en ë "u LÁ| LÁ 
array D „ " n en array R " Ĳ 


Direct op elkaar volgende COMMON elementen volgen in 


het geheugen ook direct op elkaar. Dus het bovenstaande 


A B 3 D(1) D(2) 
P A C R(1) R(2) 


Men kan een aantal grootheden in COMMON ook een groeps= 


voorbeeld 





naam geven (labelled COMMON). B.V. in het hoofdprogramma 


COMMON/NAAM1/E,G,H/NAAM2/A,B,C,D(2)//S,T 
terwijl b.v. een subprogramma (1) 
COMMON/NAAM2/P,A,C‚R(2) 
in een ander subprogramma (2) 
COMMON/NAAM1/U ,V ‚W 
en in een derde b.v. 
COMMON/NAAM2/A,B,C,D(2)//X,Y 


gebruikt worden. 
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Dit kunnen we weergeven in het volgend schema 


NAAM1 NAAM? blank 





E IG |H SIT 

U {VW XY 
verbinding hoofdprogr. hoofdprogramma hoofdprogramma 
tussen subprogr. 2 subprogramma 1 | subprogramma 3 


subprogramma 3 


Opmerking: Identifiers die niet tot de parameters behoren en 
ook niet in COMMON staan zijn locaal. Men kan COMMON dus ook 
gebruiken om geheugenruimte te besparen. | 


Formele parameters mogen niet in COMMON voorkomen. 
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Hoofdstuk VII. De programma-eenheid; Comment 


De z.g. niet uitvoerbare statements, waartoe behoren: 
de type declaraties (zoals: INTEGER, REAL, etc.) en verder _ 
DIMENSION, EQUIVALENCE, COMMON en EXTERNAL, moeten aan het be- 
gin staan in iedere programma-eenheid. (Onderling is geen | 
volgorde voorgeschreven). Daarna volgen dan eventueel de 
statement functies en vervolgens de rest van de programma- 
body. | 
Op de Cyber moet het hoofdprogramma beginnen met een kaart 
waarop de programmanaam gegeven wordt en waarop bovendien de 
in- en uitvoer media gespecificeerd worden. Voor het geval 
men alleen data van kaarten leest en de output naar de regel- 


drukker moet, kan men volstaan met 
PROGRAM naam( INPUT ,OUTPUT) 


naam = door U te kiezen naam van het programma (max. 7 symbolen) 


Een programma-eenheid heeft een der vormen 


a) D (Declaraties en andere niet uitvoerbare statements) 
SF (Statement Functies) 
S _ (Statements) 
END 

b) t FUNCTION f(a, ‚...,a,) 


(Declaraties, specificaties en andere niet uitvoer- 
bare statements) 

SF 

5 

RETURN 

END 


c) __SUBROUTINE s(a,,...,a,) of s 
D (Declaraties, specificaties en andere niet uitvoer- 
bare statements) 
DE 
5 
RETURN 
END 


Daar in Fortran iedere programma-eenheid volkomen op 
zichzelf staat, kan zij apart vertaald worden. Dit houdt in 
dat men van programma-eenheden waarin men niets verandert een 
z.g. binair deck kan maken, zodat bij volgende runs deze 
programma-eenheden niet meer vertaald hoeven te worden. Dit 
kan natuurlijk ook voor een geheel programma dat men steeds 


met nieuwe data wil gaan draaien. (zie ACCU MAP nr. 3). 


VII.2 De commentkaart (Ut) 


Een commentkaart is een kaart waarin in kolom 1 een C 
is geponst. Ze dient om commentaar te geven in de programma- 
tekst en wordt buiten beschouwing gelaten bij de compilatie 
van het programma. Het comment mag beginnen in kolom 2 en 


doorlopen tot en met kolom 80. 
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Hoofdstuk VIII. In- en uitvoerstatements (1/0) (11) 


In tegenstelling met Algol 60 behoren de 1/0 statements 
in Fortran wel tot de officiële standaard. Dit is gedaan om « 
de uitwisselbaarheid van programma's te vergemakkelijken. De 
in- en uitvoer in Fortran verschilt echter nogal veel van de 
in- en uitvoer die in Algol gebruikelijk is. De ervaring heeft 
dan ook geleerd dat dit hoofdstuk de Algol-kenners meestal de 
meeste moeite kost. We zullen daarom zeer eenvoudig beginnen 
en langzamerhand wat ingewikkelder 1/0 mogelijkheden gaan 
behandelen. | | 


De statement om invoer te plegen is, 
READ f‚k 


met f de label van de statement waarin het formaat van het 
input-veld wordt beschreven. 
k de lijst van variabelen, waaraan de input-waarden moeten 


worden toegekend. 


Als we b.v. voor de integer variabele J een waarde willen in- 


lezen, kan dat gebeuren met 


READ 1,J 
1 _FORMAT(110) 


READ 1 heeft tot gevolg dat een nieuwe kaart wordt gelezen 
volgens de specificaties van de FORMAT-statement met label 1. 
De READ is dus een uitvoerbare statement terwijl de FORMAT- 
statement een niet uitvoerbare statement is, die op elke 
willekeurige plaats in de betreffende programma-eenheid mag 
Staan, In een datd-kaart mogen älte 8) kolommen data bevatten. 

De "format descriptor!" I10, die het input-veld beschrijft, 
geeft aan dat het eerste getal geponst is In de eerste 10 
kolommen en de representatie van een integer constante heeft. 
Deze integer-waarde moet nu worden: toegekend aan de variabele 
J van de lijst. 


Op analoge wijze kunnen we getallen uitvoeren. 
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De uitvoering van de statement 


PRINT 1,J 
1 FORMAT (1Hb,120) 


heeft tot gevolg dat de waarde van J geprint wordt via de regel- 
drukker op de eerste 20 posities van een nieuwe regel. Dit 
format begint met 1Hb, waarbij we met b aangeven dat de kolom 
het symbool blank (=spatie = geen ponsing) bevat. Deze extra 
toevoeging in vergelijking met het inlezen is nodig bij het 
afdrukken via de regeldrukker, omdat dan het eerste symbool 


van iedere regel gebruikt wordt om de regeldrukker te besturen. 


„Met 1Hb zorgen we ervoor dat aan het begin van een nieuwe 


regel gestart wordt met het printen. In paragraaf VIII.9 komt 
dit onderwerp nader aan de orde. We willen hier alleen op- 
merken dat we dus nooit hetzelfde format zullen gebruiken om 
te lezen en te printen en dat we in een format dat het printen 
beschrijft steeds aan het begin van een regel een 1Hb zullen 


plaatsen. 
Voorbeeld: 


Een programma moet 100 getallen lezen. De getallen staan Op 
kaarten. Op elke kaart één getal (integer) en wel in de kolom- 
men’ 1 tot en met 10. Het programma moet het kleinste en groot- 


ste van deze getallen bepalen 


I=1 
READ 2,M 
2 _FORMAT(I10) 
MKL=M 
MGR=M 
10 READ 2,M 
IF(M.LT.MKL) MKL=M 
IF(M.GT.MGR) MGR=M 
I=I+j 
IF(I-100) 10,11,11 
11 PRINT 3,MKL 
PRINT 3,MGR 
3 FORMAT (1Hb,120) 
END | 


ne nde atd a: 
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In dit geval worden MKL en MGR onder elkaar afgedrukt. 
Willen we ze op dezelfde regel hebben dan kunnen we het 


laatste gedeelte van het programma vervangen door 


11 PRINT 3,MKL,MGR 
3 _FORMAT(1Hb,2120) 
END 


Hierbij komt MKL dan in de eerste 20 posities en MGR in 


de tweede 20 posities van de regel. 


Beschouwen nu nog eens wat nader het inlezen van integer 


getallen in het bovenstaande voorbeeld; 


| | | waarde M 
inputkaarten kolom | na READ 
1 2 3 U 5 6 7 8 9 10 11 
le kaart b hb b kh b b b b 3 +3 
2e kaart b B hb hb be bb A bh B +300 
3e kaart b == bbh B B be hb 4 OQ -30 
be kaart b bb + bb 4 0 bb +4 000 
5e kaart b bb B bk Bb hk Db B k PD 0 


Zoals al eerder opgemerkt, geven we met b aan dat de kolom het 


symbool blank (= spatie = geen ponsing) bevat. In numerieke 


data worden blanks als nullen opgevat. 


Staat de informatie niet in kolom 1 t/m 10 maar b.v. in kolom 
31 t/m bO dan kunnen de symbolen van 1 t/m 30 geskipt worden 
door de FORMAT-statement 


2 FORMATL3OX, 140) 


De descriptor nX zal bij input, n symbolen skippen. Bij output 


zullen n spaties worden toegevoegd. (n moet een positieve 
integer constante zijn). 
We hebben al gezien dat de lijst uit meerdere variabelen 


kan bestaan b.v. 
READ 10OsLediskslsMN 


Met deze statement willen we 6 getallen lezen en toekennen aan 


_ 3Ì =— 


resp. I, Js, K, L, M en N. Voor elk lijst-element is in de 


FORMAT-statement één descriptor nodig. Met 
10 FORMAT(I6,I5,213,12,11) 


is de volgende relatie vastgelegd: 


lijstelement | format descriptor 
16 
bis, 
13 
13 
EZ 
Ij 


AS A HH 


Hieruit blijkt dus dat 213 overeenkomt met 13, 13. 
In het algemeen: nIw is n maal Iw. (n positieve integer 


constante) 


Bevat de te lezen kaart de volgende symbolen, 


4 23 Hb 5 6 7 & 9 10 1 12 13 AR -L5 16 17 14-19 20 21 





dan vindt het volgende transport plaats: 


naar — waarde van lijst- 


kolommen symbolen lijt-element element na READ 
1 t/m 6 123016 E +123016 
7 t/m 11 50123 J +50123 

12 t/m 14 456 K +456 

15 Elm 17 132 L £ LZ 

18 t/m 19 12 M +12 

20 6 N +6 


Verandert men de FORMAT-statement In 
10 FORMAT(6I3) 


dan worden aan I, J, K‚, L, M en N resp. de waarden 
+123, +16, +501, +234, +561 en +321 toegekend. 
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We hebben hiervan al gezien dat voor het inlezen van 
integer data de veld-descriptor Iw gebruikt wordt, waarbij w 
de veld-breedte aangeeft (w moet een natuurlijk getal zijn). 
Het overeenkomstige lijstelement moet van het type integer 
zijn. De betreffende symbolenrij (data) moet de vorm van een 
integer constante hebben. 


Fout is dus? 


READ 10,IC 
10 FORMAT(IS5) 


met data-kaart 


want hierin heeft 30.3 namelijk door de punt in kolom 4 niet 


de vorm van een integer constante. 


Bij het printen via de regeldrukker zal de uitvoering van de 


statement 


PRINT 10.6 
10 FORMAT(1Hb,I8) 

tot gevolg hebben dat de waarde van IC wordt afgedrukt. Als 
IC b.v. de waarde +37 heeft, wordt dit bbbbbb37. Als het aantal 
significante cijfers kleiner is dan w, worden er aan de linker- 
kant spaties toegevoegd. Voor negatieve waarden komt het min- 
teken direct voor het eerste significante cijfer. Dus als IC 
de waarde -125 heeft, wordt er volgens bovenstaande format 
bbbb-125 geprint. 

Het aantal uit te voeren symbolen mag niet groter dan w 
zijn. Is dit toch het geval dan worden op de Cyber w-1 laatste 


symbolen afgedrukt voorafgegaan door een «#. 


Voorbeelden: 


PRINT 10,J 
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waarde J descriptor output 
+12345 | 18 bbb12345 
-12345 EE: b-12345 
+12345 16 b12345 
-12345 I5 “23u45 
+12345 I5 12345 


Deze descriptoren regelen de in- en uitvoer van lijst- 
elementen van type real, double precision en complex. De E, 
F en G descriptoren hebben betrekking op elementen van het 
type real of complex, waarbij voor een complex getal twee 
descriptoren nodig zijn; Één voor het reële deel en één voor 
het imaginaire deel. De D descriptor moet gebruikt worden 
voor lijstelementen van het type double precision. 


De vorm van de descriptoren is: 
Ew.d Fw.d Gw.d en Dw.d 


w: de veldbreedte, geeft het aantal externe symbolen aans; 
w moet een natuurlijk getal zijn. 
d: is een natuurlijk getal, dat het aantal cijfers achter 


de decimale punt in de externe representatie aangeeft. 


Bij input hebben de E‚ F en G descriptoren dezelfde betekenis, 


maar bij output geven ze verschillende resultaten. 


Input 


Bij input bestaat de meest algemene vorm van de externe 


rij symbolen uit: 


1) voorafgaande blanks 

2) teken + of - (geen teken is +) 

3) een rij cijfers basis vorm 
hj) een punt 

5) een rij cijfers 


en di ee een 


Eeen een 








ee 5 


6) het symbool E (of D voor double precision grootheden) 
7) teken + of = (blank of geen teken is +) 


8) rij cijfers 


Elk van deze punten kan naar keuze worden weggelaten, 


maar men moet er wel rekening mee houden dat er altijd w 


kolommen gelezen worden. Is de punt niet aanwezig in de basis=- 


vorm, dan geeft de d in de descriptor aan dat de laatste d 
cijfers van de punten 3 plus 5 achter de decimale punt ge- 
plaatst zullen worden. M.a.w. de punt komt te staan voor de 
laatste d cijfers van de basis-vorm. Staat er wel een punt 


in de basis vorm dan wordt de d specificatie genegeerd. 
Voorbeeld 1: 


READ 10,A 
10 FORMAT(F5.1) 


Data=-kaart 


A krijgt de waarde 12.34. De punt op de datakaart over- 


heerst. 


Voorbeeld 2: 


READ 10,A,P 
10 FORMAT(F5.1,E10.3) 


Datakaart 


Ll 2 4 B 5 B A 8 A 10 AL 12 AA LU 15 





A krijgt de waarde 1234.5 
en B de waarde 123.44 





ce AÔ 


Voorbeeld 3: 
Op de inputkaart staan de volgende symbolen: 


1 2 dà H B 6 7 B 4 10 14 12 43 14 15 16 17 28 19 20 Z1 


/ 
b.3 0 z 3 3 + 3 a 0 3 2 B B B 2 … Dl 4 A E 





22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 UO U1 42 43 UU 





b 0 3 b b 3 0 3 . 3 3 D - 2 bb bb bE + 2 b 


Uitvoering van de statement 
READ 10,A,B,C,D,E | 
10 FORMAT(F6.4,E7.0,G11.4,D11.4,E8.1) 


heeft tot gevolg, dat de symbolen van de inputkaart gelezen 


worden en het volgende transport plaatsvindt: 


naar waarde van lijst- 
kolommen symbolen | lijst-element element na READ 
tim 6 b30.33 À +30.33 
7 t/m 13 +3.033b B +3: 033 
ik t/m 24 bb3.033Eb03 C +303 
25 t/m 35 bb303.33D-2 D +3,0333 
36 t/m 43 bbbbbE+2 E +0. 


Hierbij moeten A, B, C en E variabelen van type real zijn en D 
een variabele van type double precision. | 
In de eerste vier gevallen komt er een punt in de basis vorm 
voor. Deze overheerst de aanduiding van de punt in de descrip- 
tor. In het laatste geval bestaat de basis vorm uit loter 
blanks. Deze worden als nul geïnterpreteerd, zodat de waarde van 
E na de READ ook nul is. n 

Omdat in bovenstaande gevallen de decimale punt dus steeds 
aanwezig is in de basis vorm, of de basis vorm uit louter 
blanks (nullen) bestaat, zou met de FORMAT-statement 


10 FORMAT(F6.1,F7.5,G11.3,D11.0,E8.5) 


precies hetzelfde worden bereikt. 
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Nog een voorbeeld: 


Op de inputkaart staan de volgende symbolen: 


1 24 4% 5 6 7 8 9 20 14 12 23 48 15 16 17 28 19 ZO 





Uitvoering van de statement 


READ 10,A,B,C 
10 FORMAT(F7.3,F8.4,FU.3) 


heeft tot gevolg dat het volgende transport benevens toekenning 


plaatsvindt 
| naar waarde lijst- 
kolommen symbolen _jst-element element na READ 
1 t/m 7 b3456bb A +345,6 
8 t/m 15 J23HEEILD B | +1234 000 000.0 


16 t/m 19 bbbb 5 0.0 


Merk op dat het eerste getal groot +345600 door de d specifi- 
catie (van Fw.d) een punt krijgt 3 plaatsen van de rechter- 
kant af, in de basis vorm. Het tweede getal, waarvan de basis 
vorm +1234 is, krijgt door de d specificatie een basis vorm 
+0.1234. Door de exponent +1b, let wel = +10, wordt de gelezen 
waarde 1 234 000 000. Het derde veld bestaande uit louter 


blanks wordt als nul gelezen. 


Tot slot van dit gedeelte over de input in samenhang met de 
E, F‚, G en D deseriptoren volgen hier nog een aantal voorbeelden: 


waarde lijst- 


format descr. Input-symbolen element na READ 
Er d 1234b +123. 
E10.2 bbbbbbbbbb = 
El2t bbbbbbbbbbb3 é +0.0003 
B1Z.4 bbbbbbbbbb3. +3 .0 
Eda 1.00bbbbbbbO int 











- 37 = 
waarde lijst- 

format descr. input-symbo len element na READ 

E12.4 bbbbbb7bbbbb „10.0 

E12.4 bbbb2b3bb.bb +20300.0 

E12.4 bbbb1.234bE3 __ _+1234.0 

E12.4 -12345678E-3 -1.2345678 

E12.0 bbbbbb+123456 +123456. 

E7.0 bbbb1LE1 F10,0 

E10.3 bbbbbbE+06 +00 
Output 


Bij output moet de interne waarde van een Fortran groot- 
heid getransformeerd worden naar een aantal (w) externe sym- 
bolen.In tegenstelling met input zijn er nu wel verschillen 


tussen de diverse descriptoren. 
Fw.d (vergelijkbaar met FIXT in Algol) 


Hierbij worden er w symbolen geproduceerd die van links naar 


rechts bestaan uit: 


1) symbolen blank 

2) — teken indien de grootheid negatief is 
3) rij cijfers (gehele deel) 

U) punt 


5) rij van cijfers (breuk-gedeelte) 


De symbolen blank worden alleen geproduceerd als het aantal 
symbolen van de punten 2 t/m 5 kleiner dan de veldbreedte w 
is. Het breuk-gedeelte bestaat uit d cijfers. Daar de punt 
ook een positie innnemt, zijn er voor positieve grootheden | 
nog w-(d+1) plaatsen beschikbaar voor het gehele deel en voor 
negatieve grootheden (waar het minteken ook nog een plaats 
opeist) maar w-(d+2) plaatsen. Er moet dus gelden 


ow (d+2) w-(d+1) 


_1 < grootheid < 10 
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Voorbeelden: 


geproduceerde 

format descr. waarde grootheid symbolen-rij 

F6.3 +0.00123 b0.001 

F6.3 +0.12468 | | b0.125 

F6.3 +1.23456 D1.235 

F6.3 -6.00000 | -6.000 

F6.3 +123.456  _#3.456 

F6.0 +123.456 bb123. 

F6.0 +0.0 bbbbO. 
In het 59° geval is de veldbreedte te klein. Het gehele deel 


6- (341) 400. De + aan het begin 


123 is n.l. niet kleiner dan 10 
van de rij symbolen geeft op de Cyber aan dat de veldbreedte 
niet groot genoeg was. Bij de laatste twee gevallen is de 
output gegeven, die geproduceerd wordt met FIN. Onder RUN zal 
in die gevallen de decimale punt niet geprint worden en links 


een extra blank toegevoegd worden. 


Ew.d (vergelijkbaar met FLOT in Algol) 


Hierbij worden w symbolen geproduceerd die van links naar 


rechts bestaan uit: 


1) symbolen blank 


2) - teken, indien de grootheid negatief is 
4 punt 
U) rij van d cijfers 


5} symbool E 
6) teken + of -— 


7) twee cijfers 


De symbolen blank worden alleen geproduceerd als het aantal 
symbolen van de punten 2 t/m 7 kleiner dan w is. Aangezien de 
exponent (5 t/m 7) 4 posities ‘inneemt en de punt 1 positie 
moet voor positieve grootheden w>d+5 en voor negatieve 

w>d+6 gelden. Indien hieraan niet voldaan is, wordt als meest 


linker symbool een «* geprint. 


Als de exponent in absolute waarde groter is dan 99, 


krijgt het exponent-gedeelte de vorm: 


5) teken + of -— 


6) drie cijfers 


Bovenstaande beschrijving heeft betrekking op de FTN- 
compiler. Indien het programma onder RUN verwerkt wordt, moet 
voor 3) gelezen worden: meest significante cijfer gevolgd 

door een punt. | 
Dus moet voor positieve grootheden w>d+6 en voor negatieve 


Ww>d+7 zijn. 


Enkele voorbeelden (ETN): 


waarde | format- geproduceerde 
lijst-element descriptor symbolen 
+12.3456 E14. 4 bbbbb.1235E+02 
-0.0012321 E13. 4 bbb-.1232E-02 
Ur T6 E12, bb-.1760E+00 
toa E9.4 1234E4+00 
012345 E94 *.123E+00 
+0.12343 E84 *.123U3E+ 
+1.x101°7 E11. bb.1000+138 
risi 2 E11. 4 bb.1000-136 


In de voorbeelden k en 5 is de veldbreedte w te klein, 
hetgeen tot gevolg heeft dat de geproduceerde rij symbolen 


begint met een %#. 


Dw .d 





Moet gebruikt worden in samenhang met double precision 
grootheden. De output is analoog aan die welke geproduceerd 
wordt met de E descriptor, alleen wordt het symbool E ver- 


vangen door een D. 


Gw.d 


Bij ouput volgens de G descriptor is de geproduceerde 


symbolen-rij afhankelijk van de groote van het reële getal 


= HO 


dat afgedrukt moet worden. Als de decimale punt niet binnen 
de eerste significante cijfers ligt, wordt volgens Ew.d afge- 
drukt. In het andere geval wordt volgens de F descriptor 
geprint, waarbij aan het eind kk spaties worden toegevoegd 


(in plaats van E+00). 


Voorbeelden: 

waarde van geproduceerde 
lijst-element descriptor symbolen 
12.345 610,3 bb12.3bbbb 

„0, LAS 5E | G10.3 b-.123bbbb 
0.012345 | G10.3 Db. L23E- (4 
123.456 | G10.3 bb123.bbbb 
1234.56 | G10.3 bb.123E+04 


Met de schalingsfactor kan de positie van de decimale 
punt van een reêel getal beïnvloed worden. Een schalingsfactor 


mag voorafgaan aan de D, E‚ F en G descriptoren. 
nPDw.d, nPEw.d, nPFw.d, nPGw.d 


De schalingsfactor n moet een integer constante zijn. 

Indien eenmaal een schalingsfactor gespecificeerd is, 
geldt ze voor alle D, E‚ F en G descriptoren in de betreffende 
FORMAT-statement totdat een andere schalingsfactor optreedt. 

Bij input wordt het getal in het inputveld gedeeld door 


fd. 

Voorbeelden: 

inputsymbolen format descr. sch.factor waarde 
123.45 Zr b ad 2 12345 
123.45 „Ar Ded mi 12345. 
b12345 IPFB:2 id: 12,35 
bbbbi23,45 „-iPE10,3 mil 120 


bbbb123.45 SPG les 3 Od Z24ho 
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Bij output hangt het effect van de schalingsfactor weer af 
van de gebruikte desciptor. Als de F descriptor gebruikt wordt 
wordt de interne waarde met 10° vermenigvuldigd. 


bÔ 


Voorbeeld: 


De interne waarde is 3.1415926 


format descr. schalingsfactor geproduceerde symbolen 
Bl 0 _bbbbb3.142 
LPFL0,3 1 ___bbbb31.416 
SPEF10,3 3 bb3141.593 
siPEIO 3 =d bbbbbb. 314 
=3PEF10.3 _3 bbbbbb. 003 


Bij de E en D descriptoren wordt de decimale punt n 


plaatsen naar rechts (als n negatief is naar links) geschoven 
en de exponent met n verlaagd. 


Voorbeeld: 


De interne waarde is 3.1415 


format descr. schalingsfactor geproduceerde symbolen 
Eds 2 0 bb.31E+01 
iPE9.2 1 b3.14E+00 
2PE9,2 2 31.42E-01 
3PE9. 2 3 “U.,16E-02 
-1PE9. 2 mid | bb.03E+02 


Bij de G descriptor heeft de schalingsfactor geen effect 
als volgens de F descriptor geprint kan worden. Als volgens 


de E descriptor geprint moet worden geldt hetzelfde als 
daarbij beschreven. 


VIII. Descriptoren voor Hollerith.tekst (strings) 


a) Aw descriptor’ 


Met de A descriptor kan tekst gelezen of geprint worden. 





ss U 


De tekst correspondeert met de inhoud van een of meerdere 


geheugenwoorden, die worden aangegeven door de lijstelementen. 


Aangezien er geen variabelen van het type Hollerith zijn, is 
het betreffende lijstelement van een ander data type (b.v. 


integer of real). Als men teksten wil vergelijken moet men 


zorgen dat de betreffende variabelen van hetzelfde data type 


Za me 

In Één geheugenwoord kan slechts een beperkt aantal sym- 
bolen opgeborgen worden. Dit aantal zal worden aangegeven met 
g. Op de Cyber is g=10, d.w.z. dat in een array van 8 ele- 


menten een kaart met 80 kolommen ingelezen kan worden. 
Voorbeeld: 


DIMENSION IA(8) 
READ 10,IA 
10 FORMAT(8A10) 


leest een hele kaart. 


PRINT 11,IA 
11 FORMAT (1Hb,8A10) 


kan de tekst weer afgedrukt worden. Op deze wijze kan men 
variabele tekst-uitvoer in een programma plegen; b;v. nummer, 


meetserie, datum uitvoering 
Wat gebeurt er als wig is? 


Input W>g: de eerste w-g symbolen worden geskipt. de laatste 
g symbolen worden in het betreffende geheugen= 
woord opgeborgen. | 

w<g: er worden w-g symbolen blank (spaties) aan het 


eind van het geheugenwoord toegevoegd. 


Voorbeeld (g=b) 


“externe symbolen format decsr. interne symbolenrij 
Â | A1 Abbb 
ABCD | A2 ABbb_ 
ABCD Al ABCD 


ABCDEFGH A8 EFGH 





Output 


Voorbeeld 


W>g: er worden eerst w-g symbolen blank Eeproduceerd, 


gevolgd door de g symbolen van het betreffende 


lijstelement 
W<g: de meest linkse 


worden geprint. 


(g=4) 


lijstelement 


bevat de 


ABCD 
ABCD 
ABCD 
ABCD 


symbolen 


bj) ni descriptor 


Deze descriptor kan gebruikt worden om de output met 


tekst te 





symbolen van het lijstelement 


format descr. 


Al 
A3 
Alt 
Al 


symbolen output 


A 

ABC 
ABCD 
DbbABCD 


& 


larderen. Bij deze descriptor hoort geen lijstelement. 


Voorbeeld: 


LG 


PRINT 10,1 


FORMAT (1Hb,I4,7HbGULDEN) 


Als I=1234 wordt geprint 


1234 GULDEN 


De n symbolen die na de H (van OEE volgen worden 


als tekst in de output opgenomen. 


Om onverwachte effecten te vermijden is het beter om in 


een FORMAT-statement,die betrekking heeft op invoer, nooit een 


H descriptor te gebruiken. 


Op de Cyber mag de string ook tussen « tekens staan. 


maakt het precies tellen van de rij symbolen overbodig. 


het bovenstaande voorbeeld mag de FORMAT-statement dan ook 


luiden: 


10 FORMAT (#b#,IU ,#bGULDEN*) 


In 


Dat 


EET IE VOED EMO OENE EEEN SE NE TN GE OERLE OE EO OG OM 


nn ie Ek DERDEN EAN ETE EUD 


EL OE EEN ED DAENS PT NO EINE UE EDE OPE OR 


ee dek ani il 


re ee ete den RE and nnn de nn Me ae ait annie eins 
- Ei 


er dE 


De descriptor Lw wordt gebruikt bij variabelen van het 


type logical (booleans). 


Input 


Door het symbool T of F in het inputveld kan de waarde 
true of false worden toegekend aan het lijstelement. Als 
W>1 kunnen er symbolen blank voor de T of F staan en/of 


andere symbolen erna. 
Voorbeeld: 


LOGICAL A 
READ 10,Â 
10 FORMAT(L5) 


inputkaart interne waarde 
bbTHb „TRUE. 
bFDAC „FALSE. 
bTRUE TRUE, 
Output 


De geproduceerde symbolen-rij bestaat uit w-1 symbolen 
blank, gevolgd door een T of een F,‚ naar gelang het lijst- 


element de waarde true of false heeft. 


nX Bij input worden n symbolen van het input-record over- 
geslagen. (n moet een constante zijn) 
Bij output worden n spaties aan het output-record toe- 


gevoegd. 


5 
/ 
| 
| 
| 
4 





Een FORMAT-statement bestaat uit een aantal descriptoren. 
In het voorgaande waren deze descriptoren altijd van elkeer 
gescheiden door een komma. In plaats van de komma kan ook het 
deelteken (/) gebruikt worden. Maar dit scheidingsteken heeft 
een extra effect, want telkens wanneer een / in de FORMAT 
voorkomt, zal een nieuw record (een record is een kaart of 


een regel) worden gevraagd of gecreëerd. Zo zullen met 


READ 10,K,L 
10 FORMAT(I10/110) 


twee kaarten worden gelezen. 

Men mag meerdere deeltekens achter elkaar plaatsen. Elk 
deelteken vraagt of levert een nieuw record. (Denk erom 
dat bij printen ieder record moet beginnen met 1Hb). 

Zoals ook al even ter sprake is gekomen, kan men voor 
een descriptor een natuurlijk getal plaatsen, dat aangeeft 
hoe vaak de betreffende descriptor genomen moet worden. Ook 
voor een groep (dat zijn een of meerdere descriptoren en/of 
scheidingstekens, die tussen haakjes geplaatst zijn) kan 


van een zodanige teller voorzien worden. 


Voorbeeld: 


KEAD TĲ,A,BGsbDsbst 
LO FORMATI2EF5 3, 2CELO. U FPB28)) 


lijstelement toegevoegde descriptor 
A | F5.3 
B | E53 
G E10 
D EBS 
5 F10.4 
Ë Fog 


Er wordt slechts één record gelezen. 


Zou de lijst nog een element G bevatten dan zou de FORMAT 
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uitgeput zijn. Dit zou in dit geval tot gevolg hebben dat een 
nieuw record wordt gelezen en G volgens de descriptor F10.4 
gelezen zou worden. 

Als een FORMAT uitgeput is, terwijl er nog lijstelementen 
verwerkt moeten worden, wordt overgegaan op een nieuw record 
en de FORMAT wordt opnieuw begonnen bij het openingshaakje 
dat hoort bij het laatste groepssluithaakje voor het eind- 
haakje. Indien er geen groep is wordt de beschrijving weer 


opgenomen aan het begin van de FORMAT. 


De lijst van variabelen mag de namen van eenvoudige 


variabelen, arrayelementen of arraynamen bevatten. Als de 


naam van een array zonder indices voorkomt wordt het: gehele 


array gelezen of afgedrukt. Hierbij wordt het array afgewerkt 
in de volgorde waarin het in het geheugen is opgeslagen. Dat 
wil dus zeggen dat b.v. een tweedimensionaal array koloms- 


gewijze wordt behandeld. 


Voorbeeld 1: 


DIMENSION A(2,2) 
READ 10,A 
10 FORMAT(UFU.1) 


met inputkaart 


L 2 8 B 5 6 7 B, 9 10 TA 12 43 1 f5 1B 17 18 48 20 





heeft na de uitvoering van de READ statement de volgende toe- 


kenning tot gevolg: 


arrayelement. . toegekende waarde 
AC1,1) d.i 
A(2,1) | | Ball 
ALA 2) eel 
A(2,2) ded 


Er 


Wil men een tweedimensionaal array rij voor rij inlezen 


dan zal men dus anders te werk moeten gaan. B.v. 


DO 11 I=1,2 

bo 41 Jets? 
11 READ 10,A(I,J) 
10 FORMAT(FU.1) 


Het array wordt nu wel in de volgorde A(1,1), A(1,2); 
A(2,1), A(2,2) ingelezen, maar omdat voor ieder element een 
READ gepleegd wordt, worden er 4 data-kaarten gelezen. Staan 
alle input data in één record, dan mag slechts eenmaal een 
READ-statement worden uitgevoerd. Men zou dan eigenlijk de 
DO statement achter de READ, dus in de lijst willen hebben. 
Hiertoe is de "implied DO" lijst ontworpen. Dit is een verkorte 
schrijfwijze voor een DO statement in een 1/0 lijst. Het ge- 
bruik hiervan zal aan de hand van enkele voorbeelden toegelicht 


worden. 
Voorbeeld 2: 


Stel we hebben een eendimensionaal array dat we tot een 


bepaalde index met ingelezen waarden willen vullen. 


DIMENSION B(15) 
N=6 
READ 10,(B(I),I=1,N) 
10 FORMAT(10F8.3) 
READ 10,(B(I),I=7,15,2) 


Door de eerste READ worden in B(1) t/m B(N=6) waarden 
ingelezen van de eerste data-kaart. (Als N groter dan 10 zou 
zijn raakt de format-beschrijving uitgeput en wordt een tweede 
data-kaart genomen). Men zij er nogmaals op gewezen dat 
FORMAT(F8.3) in dat geval tot gevolg zou hebben dat voor ieder 
array-=element een nieuwe data-kaart gelezen zou worden. 

De tweede READ leest in B(7),:-B(9), B(11), B(13) en B(15) 


waarden van de volgende data-kaart. 





aan AIEE an 


Voorbeeld 3: 


Het tweedimensionale array A uit voorbeeld 1 wordt nu rij 


voor rij ingelezen waarbij de data op een kaart staan. 


DIMENSION 
READ 10,((A(I,J), J=1,2),1=1,2) 
10 FORMAT(UFU.1) 


met data-kaart” 


ll 2 4 4 5 B 7 8 ZG 1D 14 12 13 LH 15 16 17 





heeft dat de volgende toekenning tot gevolg: 


array-element toegekende waarde 
A(1,1) Bd: 
ACL,2) Sed 
A(2,1) U.1 
AC2,2) ene 


Bij het uitschrijven wordt het eerste karakter van een 
groep data die op een rege moeten komen niet afgedrukt; maar 
deze dient om de regeldrukker te besturen. 


Enkele belangrijke karakters zijn: 


le karakter betekenis 
blank | begin op nieuwe regel 

O sla een regel over en begin 

| op nieuwe regel 
1 | begin op 1e regel nieuwe blad- 
| zijde 
+ ga terug naar begin van dezelfde 
regel | | 


Als men dus alleen maar op een nieuwe regel wil beginnen doet 
men er verstandig aan te zorgen, dat het 1e karakter altijd 


een blank is. 





Dit kan gebeuren door 


1) 1H , aan het begin te zetten (Op de Cyber dus ook # «,) 
2) 1X ‚ aan het begin te zetten | 
3) de veldgrootte van het eerste veld zo groot te kiezen, 


dat er leidende blanks voorkomen. 


Voorbeeld: 


DIMENSION A(200),B(200) 
PRINT 3,CT,ACE),BEIJ,I=110,150) 
3 FORMATCIJ,2F10, 7) 


Dit heeft tot gevolg dat er op iedere bladzijde slechts éen 


regel komt, die achtereenvolgens begint met 10, 11 , etc. 
Een beter FORMAT zou zijn 
3 FORMAT(1H ,13,2F10.7) 
Of 
3 FORMAT(1X,13,2F10.7) 
of 


3 FORMAT(IU,2F10.7) 


- De afwerking van de in- en uitvoeropdrachten 
READ f‚k en _ PRINT f‚,k 


geschiedt volgens de beschrijving van de FORMAT-statement met 
label ff. | 


-— De besturing van de uitvoering wordt overgedragen aan dit 


FORMAT-statement en wel bij het openingshaakje. 


- Er wordt een nieuw record gelezen, of er wordt begonnen met 


een nieuw record af te drukken. 


- De formatbeschrijving wordt van links naar rechts afgewerkt. 


N 
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Voor ieder lijstelement is een I, F,‚ E‚ G, D, A of L descrip- 
tor nodig. Voor een complex lijstelement worden 2 descrip=- 
toren verbruikt. De lijst wordt van links naar rechts afge- 
werkt. Overdracht van informatie vindt plaats volgens de 


conversieregels beschreven door de descriptor. 


Met de H en X descriptoren corresponderen geen lijstelementen. 
Overdracht van informatie direct tussen de formatbeschrij- 


vinge en het in bewerking zijnde record. 


Is de lijst uitgeput, d.w.z. zijn er geen nieuwe lijstele- 
menten meer, dan wordt de besturing door het format beëindigd, 
d.w.z. de in bewerking zijnde record wordt geschreven, dan 
wel de niet in behandeling genomen zijnde symbolen van het 


record worden geskipt. 


Bij elk deelteken / in het format, wordt de in bewerking 


zijnde record geschreven, dan wel een nieuw record gelezen. 


Men mag niet meer symbolen aanbieden aan, resp. betrekken uit, 


een record, dan deze bevat. 


Komt de besturing bij het rechter sluithaakje van de FORMAT 
statement, dan wordt nagegaan of er nog een nieuw lijst- 
element is. Is de lijst uitgeput, dan vindt beëindiging 

van de format-sturing plaats. Is de lijst nog niet uitgeput, 
dan gaat de besturing terug naar het begin van het format, 


resp. het begin van de laatste groep. 


Bij printen via de regeldrukker wordt het eerste symbool 


van iedere record gebruikt om de regeldrukker te besturen. 


hebben hier de statements 


READ f‚k 
PRINT f‚k 


gebruikt. 


Volgens de standaard moet men hier feitelijk 


READ (usE) kk 
en WRITE(u,f) k 


gebruiken. 


Hierin geeft u de unit aan waarmee de 1/0 opdracht uitgevoerd 
moet worden. Hiermee kan nu ook via andere media dan kaarten 
en regeldrukker data verwerkt worden. De PROGRAM-kaart voor 

de CDC komt er echter bij deze READ en WRITE wat anders | 
(ingewikkelder) uit te zien, zodat we ons hier maar tot het 
eenvoudigste geval beperkt hebben. Voor andere 1/0 mogelijk- 


heden raadplege men het FORTRAN EXTENDED Reference Manual. 
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APPENDIX A 


Standaardfuncties 


aantal symbolische type type 
definitie parameters naam par. functie 
ABS à r 
la | 4 TABS 1 on 
CABS c F 
DABS d d 
grootste gehele getal dann ai je 
< |al met teken van a e En nk 
IDINT d L 
AMOD Pr r 
a, (mod a) 2 MOD 
DMOD d 
AMAXO0 J r 
_AMAX1 r r 
max(al,a?2,...) >? MAX O0 1 d: 
MAX 1 r RE 
DMAX1 d d 
minCal 2 sns) >? idem met MAX vervangen 
door MIN | 
integer > real s FLOAT a r 
real > integer | IFIX r 
lail met teken van a? SIGN __r r 
|a1| met teken van a? 2 | ISIGN 
DSIGN 
| a, — min(a1i,a2) | 2 DIM En E 
IDIM ì 
eerste deel van dp getal 8 __SNGL d r 
reële deel REAL c 
complex J 


imaginaire deel AIMAG te 








definitie 
real > dp 


complex uit 2 reêle 


complexe congrigatie 


“log(a) 


10jogta) 


sin(a) 


cos(a) 


tg(a) 


arcsin(a) 


areCOSs 


arctg(a) 


arctg(al/a2) 


tgh(a) 


va 


s AD 


symbolische type type 
naam par. functie 


aantal 
parameters 


1 DBLE 
CMPLX r 
4 CONJG c c 


EXP r r 
1 DEXKP 
CEXP 


ALOG r | an 
1 DLOG d d 
CLOG 


ALOG10 
DLOG10 


a 5 
a 5 


SIN r r 
a DS IN d d 
CSIN 


COS r r 
1 DCOS d d 
CCOS 


1 TAN r r 


ASIN 
ACOS r Ie 


ATAN r L 
DATAN d d 


ATAN2 L EE 
DATAN2 d d 


di TANH r r 


SORT r r 
NE DSQRT d 
___CSQRT C 


Opmerking: Op de Cyber zijn, vooral bij Fortran Extended, nog 


veel meer functies aanwezig. Men raadplege b.v. FORTRAN 
EXTENDED Reference Manual, par. 1-8. 


APPENDIX B 


… Bi 


Voorbeelden 


B1. Een programma dat 100 getallen leest en het gemiddelde 


en de gemiddelde afwijking berekent en print. 


11 


12 


13 


PROGRAM B1 (INPUT ,OUTPUT) | 
VOORBEELD(9.1) UIT CURSUS ALGOL 60 


DIMENSION X(100) 
zs0 

READ 11,X 
FORMAT(16F5.2) 
PRINT 12,X 

FORMAT (1Hb,16F8.2) 
Do A T=i,100 
Z=Z+X(I) 

G=2/100 

Z=û 

Dó 2 1=1,100 
Z=Z+ABS(X(I)-G) 
Z=2/100 

PRINT 13,G,Z 
FORMAT(1Hb,2E25.14) 
END 


B2. Een programma dat n paren getallen (a, »a,) leest en voor 


î A 2 Ë 2 2 
leder paar print (3a, = a,)/(3a, + aj) en (3a, — a, )/(3a, + a 


PROGRAM B2 (INPUT ,OUTPUT) 

VOORBEELD (14.1) UIT CURSUS ALGOL 60 
READ 2‚N | 

FORMAT (110) 

DO 1 I=1,N 

READ 3,A1,A2 ON 
FORMAT (2F10.7) 

R=F(A1,A2) 


1 


) 


= BA, 


S=F(A2,A1) 
PRINT U‚R,S 

hb ___FORMAT(1Hb,2F20.13) 
END 


C DIT IS EEN FUNCTIE SUBPROGRAMMA 
FUNCTION F(X,Y) 
P=3&Xau2 
F=(P-Y)/P+Y) 
RETURN 
END 


We krijgen hetzelfde resultaat als we R=... en S=... 


vervangen door 
CALL JAN(A1,A2,R,S) 
en het functie subprogramma door 


SUBROUTINE JAN(X,Y,F1,F2) 
Pz3Xaea 2 

Q=3#Ya 2 

Elst Peli) 
F2=(Q-X)/(Q+X) 

RETURN 

END 


B3 Een voorbeeld van het rekenen met complexe grootheden 


en van het gebruik van een statement functie. 


PROGRAM B3 (INPUT OUTPUT) 

COMPLEX AC,CC,P,‚R‚ROOT1 ,ROOT2 ,WORT1 ,WORT2 
ROOT1(P,Q‚R)=(-Q+CSQRT(Q*Q-UEPER))/(24P) 
ROOT2(P,Q‚R)=(-Q-CSQRT(Q*Q-URPER))/ (24P) 
A=1 

B=2 

AC=CMPLX(A,0.0) 

Cost 2000) 

WORT1=ROOT1 (AC ,B,CC) 


a B 


WORT2=ROOT2(AC,B,CC) 
PRINT 1,WORT1,WORT2 
1 _ FORMAT(« DE WORTELS ZIJNw/2(E14.5,E12.5)) 


END 
By, PROGRAM B4( INPUT , OUTPUT) 
ë VOORBEELD: OPLOSSEN VAN DE DIFFERENTIAALVERGELIJKING 
C DY/DX=C-2#CHX#X-24Y MET BEGINWAARDE Y(0)=0 
C ER WORDT EEN RUNGE-KUTTA METHODE GEBRUIKT 
| C DE EXACTE OPLOSSING IS Y(X)=C#X#(1-X) 
| C 
DIMENSION XX(50),VAL(50) 
COMMON C 
EXTERNAL FN 
D READ 1,C,XA,YA,DX,K‚N 
Û 1 FORMAT(4F10.5,215) 
| PRINT 2,C,XA,YA,DX,K‚N 
2 _ FORMAT(3H0C=,E10.4,2X,4H XA=,E10.4,2X,4H YA=,E10.4, 
| 1 2X,4H DX=,E10.4,2X,3H K=,1I3,2X,3H N=,13/) 
CALL RK(FN,XA,YA,DX,K‚N ‚XX VAL) 
| PRINT 3,XA,YA,(XX(I) ,VAL(I) „I=1‚N) 
3 _FORMAT(16HO X EN Y(X) ZIJN:/U(F5.2,F11.7,4X)) 
| END 
C 
SUBROUTINE RK(FXY„XST,YST,H,K‚N,XN,YN) 
| C INTEGRATIE VAN EEN EERSTE ORDE DIFFERENTIAALVERGELIJKING 
| 8 MET RUNGE-KUTTA METHODE 
| C FXY GEEFT DY/DX ALS FUNCTIE VAN (X,Y) 
| 8 XST=BEGINWAARDE VAN X 
| C YST=BEGINWAARDE VAN Y 
| C H=STAPGROOTTE 
© K=AANTAL STAPPEN TUSSEN XN(I) EN XN(I+1) 
| C N=AANTAL XN(I) 8 
8 EINDWAARDE X IS X=XST+KeNeH 
© 


| DIMENSION XN(1),YN(1) 
H2+H/ 2 








z Bi = 


X=XST 
Y=YST 

DO 2 I=1,N 

DO 1 J=1,K 
S1=HeFXY(X,Y) 
SY=H&FXY(X+H2,Y+51/2) 
S3=H&FXY(X+H2,Y+S52/2) 
Sh=H&FXY(X+H,Y+53) 
X=X+H | 
Y=zY+(S1+2#S2+2#S3+SU)/6 
XN(I) =X 

YN(ID=Y 

RETURN 

END 


FUNCTION FN(X,Y) 
COMMON CS 
FN=CS-2#CSuXan2-2aeY 
RETURN 

END 


